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Abstract 

Different families of states, which are solutions of the time-dependent 
free Schrodinger equation, are imported from the harmonic oscillator 
using the Quantum Arnold Transformation introduced in pQ. Among 
them, infinite series of states are given that are normalizable, expand 
the whole space of solutions, are spatially multi-localized and are eigen- 
states of a suitably defined number operator. Associated with these 
states new sets of coherent and squeezed states for the free particle are 
defined representing traveling, squeezed, multi-localized wave packets. 
These states are also constructed in higher dimensions, leading to the 
quantum mechanical version of the Hermite-Gauss and Laguerre-Gauss 
states of paraxial wave optics. Some applications of these new fami- 
lies of states and procedures to experimentally realize and manipulate 
them are outlined. 

1 Introduction 

In the context of the quantum free particle, the eigenstates of the Hamilto- 
nian, which are also eigenstates of the momentum operator, are not normal- 
izable. These states, the plane waves, are fully delocalized. However, it is 
customary to expand any normalizable solution of the free, time-dependent, 
Schrodinger equation in terms of plane waves, using the Fourier transform, 
building in this way wave packets which represent localized solutions and 
that are no longer eigenstates of the Hamiltonian. The simplest examples 
are the Gaussian wave packets, which have the property of minimizing the 
uncertainty relations between the position and the momentum operator. 

In pQ a transformation was proposed that maps states and operators 
from certain quantum systems to the free particle. This transformation is 
the quantum version of the Arnold transform (QAT), which in its original 
classical version [2] maps solutions of a certain type of classical equation of 



motion, a nonhomogeneous linear second order ordinary differential equation 
(LSODE), to solutions of the classical equation for the free particle. 

In this paper we construct different families of states in the space of solu- 
tions of the quantum free particle which are the image under the QAT of the 
corresponding families in the quantum harmonic oscillator. Among them, 
we provide discrete bases of free particle states mapped from the eigenstates 
of the harmonic oscillator. They are not eigenstates of the free particle 
Hamiltonian, i.e. they are not stationary states, but, rather, eigenfunctions 
of a certain quadratic operator N with discrete eigenvalues. The first states 
of these bases are Gaussian wave packets, and are localized with arbitrary 
initial size, related to the oscillator frequencies and chosen at will, and have 
initial minimal uncertainties. The following ones are "multi-localized" in the 
sense that, for instance, in one dimension the n th -order state presents n ze- 
ros and n + 1 humps which spread out with time. This mimics the situation 
for the harmonic oscillator to such an extent that, in one dimension, it is 
possible to build creation and annihilation operators a) and a. The number 
operator proves to be N ~ a) a, although in the case of the free particle it is 
not the Hamiltonian. Going even further, we give a set of "coherent" and 
"squeezed" states which are interpreted as traveling, squeezed wave packets. 
This construction can easily be generalized to higher dimensions in different 
coordinate systems. 

These families of states and their construction through the QAT can 
be of physical relevance in the description and preparation of initial states, 
and even time-evolution processes, in Atom Optics, ion traps physics and 
Bose-Einstein Condensates. For instance, the process of turning off and on 
a harmonic potential trapping a set of atoms, as well as the production of 
coherent and squeezed states in this context, are easily achieved by means 
of the QAT. It might also be useful in describing scattering processes in a 
discrete basis, instead of the usual approach employing plane waves. 

The content of the paper is as follows. Section [2] is devoted to the con- 
struction of the discrete bases of states in one dimension by means of the 
Quantum Arnold transformation, mapping the eigenstates of the harmonic 
oscillator to the free particle Hilbert space. Section [3] deals with the con- 
struction of coherent and squeezed states and the interpretation of them as 
traveling, squeezed wave packets. Section [J] generalizes this construction to 2 
and 3 dimensions in Cartesian, polar and spherical coordinates, giving rise to 
Hermite-Gauss, Laguerre-Gauss and spherical-Gauss wave packets, respec- 
tively. Section [5] presents some possible physical applications, and proposes 
experimental settings for producing and manipulating these states. Sec- 
tion O gives some possible generalization of these constructions and outline 
possible lines of research. Finally, an Appendix collects the main formulas 
for the computation of coherent and squeezed states for the free particle. 
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2 Discrete bases of wave packets in one dimension 



Let H be the Hilbert space of solutions of the free particle, time-dependent 
Schrodinger equation, and %ho the one corresponding to the harmonic oscil- 
lator of a given frequency oj. We shall denote by tp(x, t) £% the free particle 
solutions and by ip'(x',t') £ Hho the harmonic oscillator ones. Then, the 
QAT [1], relating solutions of time-dependent Schrodinger equations for the 
free particle and the harmonic oscillator, or rather, its inverse, is given by: 

4,'{x\t') = A'^(x,t) = -^==e* ft ^ ip( 777T' ^U) > W 

where the classical Arnold transformation [2j is explicitly 

A: RxT' — >RxT 

(*',*') (x,t)=A((x',t')) = (^^). 



ttlft'K ( 2 ) 



T and T' are open intervals of the real line containing t = and t' = 0, 
respectively, ui(t') and u 2 {t') are two independent solutions of the LSODE 
(here dots mean derivation with respect to t'): 

x' + f(t')x' + co{t') 2 x' = 0, (3) 

and W(t') is the Wronskian W(t') = iii{t')u 2 (t') - ui(t')u 2 (t') of the two 
solutions. For the case of the harmonic oscillator / = and cj(i') = cj, 
and the two independent solutions can be chosen (see pQ for details) as 
u\(t') = ~ sin(cui') and u 2 (t') = cos(wt'), with W(t') = 1. It can be checked 
that the change of variables results irfl: 

t' = — tan _1 (a;t) 

x = cos(tan~ 1 (cut))x = , . (4) 

Diagrammatically, the QAT can be represented as: 

A 

n - n uo 



HO 





1 The change of variables in the particular case of the harmonic oscillator was already 
given by Jackiw [3]. 
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where U and C/ho stand for the evolution operators of the free particle and 
the harmonic oscillator, respectively, while %q is the Hilbert space, either for 
the free particle or the harmonic oscillator, of solutions of their respective 
Schrodinger equations at t = (see pQ for the conditions under which Hq is 
common to both systems). 

Thanks to the commutativity of this diagram, and the unitarity of the 
operators appearing in it, we can map objects (wave functions, operators, 
expectation values, uncertainties) from one system to the other. In pQ we 
benefited from this fact for transporting the simplicity of the free particle to 
more involved systems, finding, for instance, an analytic expression for the 
evolution operator of the complicated systems (even with time-dependent 
Hamiltonians) in terms of the free particle evolution operator. In this paper 
we shall proceed the other way round, transporting objects and properties 
from the harmonic oscillator to the free particle. 

Applying now the QAT to the time-dependent eigenstates of the har- 
monic oscillator Hamiltonian H^q, with energy E n = hw(n + 5), 



(5) 



where N n = (7^) 4 ^ 2 "n! ' we °^* a ^ n ^ ne following set of states, solutions of 
the Schrodinger equation for the free particle: 

x 2 S* i 

*.(,,.)-* ;> e W(£p*. ( * (6) 



ISf \\S\J y V2L\S\ 

where, in order to obtain a more compact notation, we have introduced the 



quantities L = y 2mZJ' w ith dimensions of length, and r = 2m ^ = uj^ 1 , 
with dimensions of time. We also denote by 5 the complex, time dependent, 
dimensionless expression 5 = 1+iut = ^+i 2 mL i = H"^t - ^ e nave a l so used 
the fact that e~ iujt ' = e^ itan_1 ( wt ) = Note that with these definitions 

the normalization factor M n can be written as M n = / ^ ■ 

This set of states constitutes a basis for the space of solutions of the 

free Schrodinger equation, since it is mapped from a basis for the harmonic 

oscillator through A, which is unitary. 

The family of wave functions ([6]) has been known in the literature as 

Hermite-Gauss wave packets [3], and they have been widely used, in their two 

dimensional version (see Section |3J, in paraxial wave optics [5p , However, 



2 In [6] a basis of discrete states and their corresponding coherent states were con- 
structed for the damped harmonic oscillator, which, in a certain limiting process, would 
reproduce ours. Also, in [7] the harmonic oscillator wave functions are mapped to solu- 
tions of generalized, time-dependent harmonic oscillator in the framework of [8], which 
resembles the approach in .1 , although with somewhat less physical insight. 
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this kind of states and the ones constructed in Sections [3] and |4] are better 
understood in the framework of the QAT. Note that, making use of the 
classical solutions only, and through the QAT, we have been able to import 
the time evolution from the stationary states of the harmonic oscillator, 
ip' n (x',tf), to the non-stationary ones, ifj n (x,t), without solving the time- 
dependent Schrodinger equation. 

The first state of this basis, the one mapped from the harmonic oscillator 
vacuum state, is given by: 

i x 2 S* _ L x 2 

5 e"^ J =^Tle~4l?S, ( 7 ) 

which is nothing other than a Gaussian wave packet with center at the origin 
and width L. The parameter r is the dispersion time of the Gaussian wave 
packet (see, for instance, [9]). 

Figure [T] shows some of these wave functions and how they evolve in 
time. 



— '0 




Figure 1: Spreading under time evolution of wave functions tpo, ipi and ip2, 
with tk = far. 

We see that the number of "parts", or humps, of the wave functions, 
determined by the number of zeros, is "quantized", in the sense that there 
is only one hump between two consecutive zeros. They are therefore multi- 
localized, as the probability of finding the particle is "peaked" in separate 
regions of space, with nearly zero probability of finding it in between. 

The QAT also allows to map invariant operators from one Hilbert space 
(Hno) t° the other (T~L). These invariant operators are also known as con- 
stant or integral of motion operators, in the sense that their matrix elements 
are constant with respect to their corresponding time evolution, and preserve 
their respective Hilbert spaces. Particularizing the general expression given 
in [I] or by direct computation, the conserved position and momentum op- 
erator for the harmonic oscillator can be written: 

X' 
P 

5 



, . . (2vr)-4 /S*\ 



ui(t') , ih , d , , ih . . d 

— — — x H ui(t )— = cosujt x H smut — 

Wyv) m ox' mu ox' 

—%riU2\t — m———-x = —incosut — — - + mu smut x , 
ox' W(t') ox' 



and therefore, the conserved creation and annihilation operators are: 



Note that these are the operators acting on solutions of the time-dependent 
Schrodinger equation for the harmonic oscillator vp' n (x',t') as ladder opera- 
tors. 

Position and momentum operators X' and P' are mapped into opera- 
tors representing conserved position and momentum operators for the free 
particle through the QAT: 

ih d 
X = x + —t— 
m ox 

ox 

As a consequence, ladder operators for the harmonic oscillator can be 
mapped into ladder operators for the free particle that act as creation and 
annihilation operators for the (time-dependent) Hermite-Gauss states: 

r r ° X 

The action of d and d^ on the Hermite-Gauss wave functions is the usual 
one: 

aip n (x,t) = ^/nil> n -i(x,t), aJip n (x,t) = Vn + lip n+ i(x,t). (9) 

It is possible to introduce this discrete basis without resorting to the QAT 
in a very intuitive manner. The key point is that the operator a annihilates 
the Gaussian wave packet, and this fact characterizes it. The whole family of 
states (J6j) can be generated by acting with the adjoint operator d' of d. The 
rest of the construction, i.e. uncertainties, coherent and squeezed states, 
etc. would proceed without the need of resorting to the QAT. However, 
as we shall see below, the QAT is very useful when performing involved 
computations in a very easy way. 

In fact, a very useful property of the QAT is its unitarity, implying that 
it preserves scalar products and therefore expectation values. This will allow 
us to compute expectation values of operators in the free particle in terms 
of the corresponding ones in the harmonic oscillator. 
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Denoting by (•, •) the scalar product in the Hilbert space of solutions of 
the free particle Schrodinger equation, and by (-,•)' the corresponding one 
in the harmonic oscillator Hilbert space, we have that: 

^A) = ^'A')' (0)^, = (0%, (10) 

where ip, 4> and O are related to ip', 4>' and O' through the QAT, respectively. 
It is important to note pQ that (|10p is only valid when O and O' are in- 
variant operators. These operators are essentially X and P above, closing a 
Heisenberg-Weyl algebra (or X' and P' in the harmonic oscillator), and their 
powers. Among them, the most important ones are the quadratic operators, 
which together with the Heisenberg-Weyl algebra expand the Schrodinger 
algebra (see ([17]) ). 

For instance, we can apply this property to compute the uncertainties 
associated with each wave function. As a function of time, for each state 
ip n (x,i), they read: 

Ax n Ap n = (n + hh\S\ = ^\6\. (11) 



To compute the uncertainties Ax n = ^ ((x) 2 ) ^ n — {{x)^, n ) 2 and Ap n = 
V ((P) 2 )i>n ~ iip)ipn) 2 we have used the transformation of non-conserved op- 
erators when they are written in terms of conserved ones, that is: 

m mu2(t') U2\t ) 

p=P ^—^ P' = —ip' + mu 2 tx') , (12) 
\°\ 



as well as 



<(j5)V = ((f") 2 )i;=((p') 2 )i;=™wS(n + i), (13) 

where in the last equality of the equations the quantum Virial Theorem for 
the harmonic oscillator potential [8], which is homogeneous of degree two, 
has been used: \mu 2 {(x')%, n = ^((pf)^ = (H' HO Y^- 

For n = the time evolution of the uncertainty relation (jlip is the one 
which results from the usual Gaussian wave packet [S], and, among all, the 
minimal one. 

The number operator associated with the creation and annihilation oper- 
ators above will provide the position of the state in this grid of uncertainties. 
We can derive it (or map it from the number operator for the harmonic os- 
cillator) in the usual way: 



N = — I a* a + a a 



d 2 t , d 1, 



x 



2 



- |5|2L2 ^ + V ( ^ + 2 ) + 4l^ 



(14) 
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By making use of the Schrodinger equation, we can turn this operator 
into a first order one: 



N 



' - l2 d t d 1 x 2 



(15) 



the validity of which is restricted to solutions of the Schrodinger equation. 
The action of this operator is such that: 



1 



Nt/j n (x, t) = (n+ -)4> n (x, t) 



En 



1p n (x,t), 



(16) 



thus reproducing the uncertainties (up to a factor K) given in Eq. (fTTj) at 
time t = 0. 

It is quite interesting to note that this operator belongs to the alge- 
bra of the "maximal kinematical" symmetry of the free particle, i.e., the 
Schrodinger group [TO] , built out of constants of motion operators up to 
quadratic order. This symmetry is the standard Galilean symmetry (with 
generators P 2 = 2mH, P, X and the identity I) together with spatial dila- 
tions (generated by XP = \{XP + PX)) and non-relativistic "conformal" 
transformations (with generator X 2 ). These generators can be written, in 
their first-order version, as: 



X = x + ih-^-4- 

™, ox 



p2 

XP 



p 

X 

X 2 = 2i±t 2 £ + 2i±tx£- + x 2 + Hit, 
m at m ox 1 m ' 

providing the non-trivial commutation relations: 

X,P 



—2ihtQj — ifrx-r^ — y 



(17) 



X,P 2 



P,P Z 



X 2 ,P 2 



2ihP 


MhXP 



X,X Z 
P,X 2 
X 2 ,XP 



ih 




-2ihX 



2ihX 2 







(18) 


X,XP 


= ihX 


(19) 


P,XP 


= -ihP 


(20) 


P 2 ,XP 


= -2ihP 2 . 


(21) 



It is easily checked that N belongs to this Lie algebra, its relation with the 
basis above being: 



N 



-(- 

hw 2m 



P 2 H X 2 ) 



1 



Huo , 



(22) 



where H^o is the operator mapped through the QAT from the harmonic 
oscillator Hamiltonian H{iq- See [1] for the relevance of quadratic operators 
like N, but distinct from the Hamiltonian, for building bases of the Hilbert 
space (see also [TT1 151 IB1 \12\ ) . 
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3 Coherent and squeezed states 



Once we have successfully imported the eigenstates of the harmonic oscilla- 
tor, we can go further with this scheme and map different families of coherent 
states from the harmonic oscillator to the free particle (as we shall see, they 
might be defined equivalently, without resorting to the QAT, as eigenstates 
of the annihilation operator a in the free particle). Even more, squeezed 
states may be introduced in a similar fashion. These states have the prop- 
erty that, having minimal uncertainty relations, the uncertainty in position 
(or momentum) is squeezed, at the expense of increasing the uncertainty in 
momentum (or position). These states are used in every day experiment in 
Quantum Optics since they have low noise and allow to measure quantities 
with higher precision |13j . 

The most general squeezed and displaced number states |n, £, a)' for 
the harmonic oscillator are built using the squeezing operador = 
e 2(£* a ' 2 ~£( a ' ) 2 ) an d the displacement operator D'(a) = e aa ' f - a * a ' ■ 

\n,Z,a)' = D'(a)S'(0\n)' (23) 
where a is, as usual, the complex number 

frfuJ 1 x .p L I 

a = y^n xo+ l v^ po = 2T + 1 ~ = 2L {xo + woT) ' (24) 

with vq = ^,and £ is the squeezing parameter that will be taken real, £ = r, 
for simplicitjo 

The complete expression for the configuration space wave functions of 
these states for the harmonic oscillator, and its free particle version through 
the QAT, is given in the Appendix. Here we only give their general form 
for the free particle: 

r M M=X n e r / 2 -j= (jJi J ^^^(g), (25) 

where q = | and 5 r = 1 + ie uit, and the phase 9[x,t) is given in 

Eq. (0SD in the Appendix. 

Now we focus ourselves on some particular cases, giving their physical 
interpretation in the free particle system. By putting n = r = in (|25|) . we 

3 The difference between real and imaginary squeezing parameter is easily seen checking 
the action of the squeezing operator on the Wigner quasi-probability distribution, with 
domain the phase space (x,p). Real £ implies squeezing in the x or p direction, while 
imaginary £ squeezes in the diagonal x + p or x — p directions. In other words, when 
£ = pe 10 , the squeezing takes place in the direction | in the plane (x,p). 



9 



obtain coherent states for the free particle: 

, i (a; - x ) 2 + x t;oi + iTv (v t -2x + x ) 

(2-k) 4 / S* \ 2 



(26) 

which could also have been obtained as eigenstates of the annihilation op- 
erator since they verify: 

a<p a (x,t) = a<j> a (x,t) . (27) 

These states can also be obtained by the action of the imported displace- 
ment operator D(a) = e aat-Q * a = e^ p ° x ~ x ° p \ which is, up to a phase, a 
Galilean boost with parameter po together with a translation of parameter 
xq on the Gaussian wave packet. They constitute an over-complete set of 
the Hilbert space of the free particle and represent traveling Gaussian wave 
packets, with mean momentum and initial position po and xq, respectively. 
These states are not eigenstates of the number operator N. The expectation 
values of this operator on these coherent states can be computed importing 
the result from the harmonic oscillator using (|1U|) and {\22f\ : 

(0 a \N\0 a ) = h(\a\ 2 + ±). (28) 



Setting now n / 0, r = in f)25|) . we obtain a family of states of the 
free particle that could be obtained by acting with Galilean boosts and 
translations on a fixed state of the basis, ijj n (x,t), and constitutes a new 
over-complete set of states for the free particle Hilbert space: 

,n, ^ K r 1 f5*Y + 3 ,X-X -V t, 

<t>a{X,t) = M n —7= [rjr] H n ' 



^W\\\5\J " v V2L\5\ ' 

(x - Xq) 2 + X V t + iTV (v t - 2x + Xq) 

t 4L 2 5 , (29) 

representing traveling multi-localized wave packets bearing n + 1 humps, 
with mean momentum and initial position po and xq, respectively, where a 
is given by (|24|) . As in the case n = 0, these states are not eigenstates of 
the number operator, but the expectation values, computed as before, are: 

^ a \N\r a ) = h(\a\ 2 + n+ 1 -). (30) 

Being a set of coherent states, the uncertainty relations of if)™, Va G C, are 
the same as those of ip n given in Eq. (fTT]l (see [13])- 

Finally, let us discuss the simplest case of squeezed state, that with r^O, 
n = and a = 0, corresponding to the squeezed vacuum in the harmonic 
oscillator or a "squeezed" Gaussian wave packet in the free particle: 
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x 2 e 2r 5* 

This wave function corresponds to a Gaussian wave packet where the 
position and time variables have been rescaled, in addition to the wave 
function, as 

tp r (x,t) = e r /Vo(e r x,e 2r t). (32) 

Again, this can be interpreted as the action of the imported displacement 
operator S(r) = e^ a ~' a > ' = e~ which turns out to be a dilation. 

Comparing the expressions of ([7]) and (|3ip . we realize that (|3ip corre- 
sponds to a Gaussian wave packet with initial width L r = e~ T L. This can 
be easily understood since in the free particle there is no natural frequency, 
and therefore there is no natural width L, so that squeezing can be absorbed 
in the width L of the wave packet. 

4 Discrete bases of wave packets in higher dimen- 
sions 

The generalization of the previous construction to higher dimensions is im- 
mediate. In more dimensions, due to the symmetries of the harmonic oscilla- 
tor, we can choose different separation of variables to solve the Schrodinger 
equation for the harmonic oscillator, which amounts to find simultaneous 
eigenstates for the harmonic oscillator Hamiltonian and other operators com- 
muting with it and among themselves. For instance, separation of variables 
in cartessian coordinates in N-dimensions is equivalent to diagonalizing si- 
multaneously the harmonic oscillator Hamiltonian, and the ID harmonic 
oscillator in each coordinate. This can always be done even if the frequen- 
cies for each direction are different. The common basis of eigenstates is 
the product of ID harmonic oscillator eigenstates in each variable. If we 
impose spherical symmetry, then all frequencies must coincide and we can 
search for common eigenstates of the Hamiltonian and the angular momen- 
tum operators commuting with it and among themselves (for instance, in 
three dimensions they would be 1? and L z , and in two dimensions it would 
be only L). 

For the case of Cartesian coordinates, the generalization of the previous 
sections is immediate by a direct generalization of the QAT to higher di- 
mensions, the discrete basis of states for the N-dimension free particle being 
products of N copies of ([6]) in each of the N variables (x±,X2, ... ,xn). If 
' l Pn i { x i^) represents a wave packet of the form given in ([6]) but in the vari- 
able X{ and with width L, (and with the corresponding frequency cjj), then 
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Figure 2: Density plot of the probability distribution for the Hermite-Gauss 
states ip(i o) and ipnu- The positions of the maxima are shown. 



the iV-dimensional wave packets are: 

ipfi{x,t) = Jlf =1 i! ni {xi,t) , (33) 

where n = (ni, n%, . . . , njy)- If all the widths Li = L are identical, it can be 
written as: 

(34) 

Coherent states, or traveling wave packets, are defined similarly, as prod- 
ucts of ./V copies of (f2Uj) , with a vector a £ : 

^,t) = Il£ 1 ^(x i ,t). (35) 

For N = 2 these states have been widely used in the paraxial approx- 
imation to the Helmholtz equation of wave optics, known as the Hermite- 
Gauss states [5j. Due to the analogy of the paraxial approximation to the 
Helmholtz equation (with the z coordinate acting as time) and Schrodinger 
equation in two dimensions, they have also been exploited in atom optics 
and matter waves |15j . 

In Figure [2] a density plot of the probability distribution for the Hermite- 
Gauss states "0(1,0) an d ^(1,1) is shown. 

In the case of cylindrical symmetry in paraxial wave optics, the polar 
version of these estates have been used, known as Laguerre-Gauss states [5j. 

The corresponding states for the free Schrodinger equation in two dimen- 
sions, derived throught the QAT from the 2D harmonic oscillator in polar 
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coordinates, are: 

7iW L " ( 2LW ) ' (36) 

where n, / = 0, 1, 2 . . ., and L l n (x) are Laguerre polynomials. The state with 
n = 0, 1 = is the Gaussian wave packet in two dimensions. 

These states are eigenstates of the angular momentum operator L in 2 
dimensions, with values Lip^i(r,(f>,t) = Htp^ir, <f>, t). In Figure [3] density 
plots of the probability distribution for the Laguerre-Gauss states an d 
i^ii are shown. 

The generalization to three dimensions in spherical coordinates is straight- 
forward, the states having the form of (|36|) but in terms of spherical har- 
monics and confluent hypergeometric functions M(-, ■ ; •) (which are also 
polynomials in this case): 

r(; + 3/2 + n-l) /j*x2(n-l)+/+3/2 
V / 2(n-l)!r(/ + 3/2)2L3|5| \W\ J 

Yr{9 ^ (y^) ' M{ ~ n + M + 3/2 ' (37) 

None of these states are eigenstates of the free particle Hamiltonian, but 
can be seen to have expectation values of the energy equal to half the energy 
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of the corresponding, through the QAT, harmonic oscillator eigenstate: 



\E H E r -i = MEili n i + f ) 

\E nl E nl = fiw(2n + Z + l) (38) 

\E nlm E nim = hu)(2(n - 1) + I + |) 

These states, as in the one dimensional case, have many humps. For 
Cartesian coordinates they have Jlf = Ani + 1) humps, while in polar coordi- 
nates they have (n + 1) humps (seen in the radial coordinate) with annular 
form. 

Uncertainty relations can also be computed for these states. However, 
only the cartessian version of the operators x and p can be used since the 
canonical momentum associated with the radial coordinate, p r , is not self- 
adjoint. The expressions, computed as in the one dimensional case, are: 

(Ax i )^(Ap i )^ 

As in the one dimensional case, squeezed states, interpreted as rescaled 
wave packets, can also be introduced in higher dimensions. 



71,1 



= -\6\E ni , i = l,...N 

= (Ay^JApy)^ = ±-\5\E nl (39) 

= ■Z—\fi\Enlm 



5 Physical applications 

The theoretical relevance of these multi-localized, squeezed, traveling wave 
packets for the free particle, as well as their construction by means of the 
QAT, is of no doubt. From a more practical point of view, the preparation 
of this kind of discretized free states might be achieved by the use of a 
harmonic oscillator the potential of which is switched off at a given tim^f|. 
The vacuum state of this harmonic oscillator, when switched off, will provide 

the "vacuum" Gaussian wave packet with width L = y^^j, where m is the 
mass of the particle and uj the frequency of the oscillator. Note that the 
dispersion time r coincides with the inverse of the frequency of the oscillator. 
If the harmonic oscillator is in the n-th excited state, the (n + l)-hump state 
is obtained. To obtain traveling states, the initial state should be a coherent 
state 4> a (x,t) of eq. ([26]) for a one-hump traveling state or (f>™(x,t) for a 
(n+l)-hump traveling state. These coherent states can be obtained by acting 

4 This a consequence, again, of H10|l and the quantum Virial theorem for the quadratic, 
homogeneous harmonic oscillator potential. 

5 This idea was proposed in [27], and was named "Quamtum Sling". 
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with time-dependent classical forces on the harmonic oscillator according to 
Glauber [28] (see also [THE]). In fact, if the classical force is given by the 
potential V(x) = —f(t)x, and the initial state is the vacuum |0), then, if the 
force f(t) has finite duration or it is fast decaying in time, after a suitable 
time a standard coherent state la) is obtained with a = i- — f(u))*, where 

f{oo) is the Fourier component of f(t) in the frequency u of the oscillator. For 
an arbitrary force f(t) a time-dependent coherent state \a(t)) is obtained, 
see [ME]. 

It is possible to "capture" one of these traveling states at a later time t\ 
switching on a harmonic oscillator potential with an appropriate frequency 
oj\ that would "freeze" and keep it in an coherent state (without dispersion) , 
until the potential is switched off again at a time ti and the wave packet 
is released, traveling again as a free wave packet that disperses in time. 
This way, information can be stored temporally in this "oscillator traps", 
which can also be used to further manipulating them or even measuring 
the resulting state by means of adequate lasers. The frequency uj\ required 
to capture the dispersed wave packet should be fine tuned in such a way 
that the wave packet, at the time t±, matches an appropriate solution (with 
the same n and the same a) of the harmonic oscillator with frequency u)\. 
Using the relation between u and L (valid at t = 0), u = 2 mL 2 > anc ^ ^ e 
fact that the width of the wave packet scales as \S\, see eq. (fT3|) . then 
oj\ = 2m L^\8 1 \ i ' w h ere $1 = l + iw£i. This expression can be written in terms 
of distances using the mean velocity vq of the wave packet encoded in the 
complex number a. 

This process can be elegantly described by a sequence of QAT and evo- 
lution operators in the following way. Denote by A w ^ {t) the QAT from a 
harmonic oscillator of frequency u performed at time t, and by U' u the uni- 
tary time evolution operator for that harmonic oscillatofi Then the state 
after the process described above is (only the time dependence is explicit): 

i;(t) = f7(i,i 2 )^ 1 (i2,ii)^(ii,io)^ / (io) = 

= ^.faW^^,^)^^,*!)^^!)^,^.^'^) = ( 4 °) 

where i^'ito) is the starting harmonic oscillator state. This way, in the last 
expression only evolution operators for harmonic oscillators appear. 

Note that if the frequency of the harmonic oscillator is not modified, 
and the same frequency uj is used to capture the state in the harmonic 
oscillator trap, the resulting state will be a squeezed state with squeezing 
parameter r given by r = — log(|<5i|), which is negative. This can be seen 

6 With this notation, for t = to, Aui,t {to) is the identity (classical solutions tii(t') and 
U2(t') have been appropiately chosen, see [1 ), and t = ■ 
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as a feasible way of producing squeezing in trapped states, simply switching 
off-switching on the trap for a period of time t, resulting in a squeezing 
parameter r = — | log(l+o; 2 t 2 ). In fact, a similar way of producing squeezed 
states in Bose-Einstein Condensates was reported in |29j. 

The states here constructed are rather robust in the sense that the num- 
ber of humps is conserved even in the presence of small perturbations. Nu- 
merical calculations in one dimension have been performed, simulating per- 
turbations by square potentials (well or barriers), leading to the conclusion 
that this holds as long as the mean energy of the state is large compared with 
the scale of the perturbing potential and the wave packet is sharp enough 
in momentum space in such a way that the transmission coefficient can be 
considered a constant. Under these circumstances (see for instance [26J), 
the wave packet behaves as a plane wave and the effect of the barrier in 
the transmitted packet is an overall attenuation, preserving its shape, and 
a time delay which takes its maximum values for energies near the resonant 
ones (and where the transmission coefficient is one). As shown in [26J, this 
result is valid for any bounded potential of compact support, provided that 
the width of the potential is small in the sense that the time to pass through 
the barrier is smaller than the dispersion time of the wave packet r. There- 
fore, the conclusions obtained with the square potential can be generalized 
to any finite-range bounded potential. 

These wave packets are also robust under the influence of time depen- 
dent, homogeneous external fields, or even linear damping. In these cases, 
the centroid of the wave packets follows the classical trajectories, but their 
shapes are unaltered, apart from the unavoidable dispersion^ [3]. The case 
of linear damping is interesting due to the fact that the presence of damp- 
ing prevents the dispersion of the wave packets, which asymptotically have 
finite width [61 H] . The case of homogeneous external fields is particularly 
interesting, since it includes the case of free fall, which is a common exper- 
imental situation in atomic physics or Bose-Einstein Condensates (see for 
instance |30j). 

Under the conditions commented above, these wave packets evolve with- 
out distortion even in the presence of perturbations. However, one could 
be interested, acting with appropriate potentials, in obtaining transitions 
between wave packets with different number of humps, in such a way that, 
for instance, a one-hump packet splits into a two-hump packet or a two- 
hump packet coalesces into a one- hump packet. This would open the door 
to performing quantum gates acting on q-bits realized with the one-hump 
and the two-hump states. A way of implementing this is to benefit from the 
fact that wave packet dynamics is similar to wave optics in the sense that 

7 Regarding the dispersion effect, which should be taken into account for practical 
application in Atomic Physics or BEC, it takes place in a time scale given by the dispersion 
time t. For instance, for a condensate of Na atoms with size « 1mm, this time scale is 
« 500s. 
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an analogue to the ABCD law for optics is satisfied for wave packets |16j . 
Even the transmutation of Hermite-Gauss wave packets into Laguerre-Gauss 
wave packets can be achieved using ABCD matrices, implementing a mode 
converter |31j . 

Among further theoretical applications, we could think of expanding 
plane waves in terms of the discrete basis {<j>~}fi£N N i with fixed a 6 C w , 
and describing scattering process in a discrete basis, or expanding arbitrary 
wave packets in a continuum over-complete set {^} JgC iv, with fixed n G 

, which could be discretized in a lattice Z N x X N of points to perform 
numerical computations, while keeping the over-complete character. Note 
that they are over-complete for t = if the volume of the unit cell is smaller 
or equal to H N , and again by the unitary time evolution they continue to be 
over-complete for any t, see |14j . Similar constructions, like the Harmonic 
Oscillator (HO) method [20] or the Transformed Harmonic Oscillator (THO) 
method [21] have been proposed, mainly in nuclear physics, to describe the 
bound states and the continuum spectrum in a discrete basis. But there 
the construction just appears as a mathematical tool for approximating the 
solutions, with no physical meaning. Our states, however, are physically 
meaningful (as traveling wave packets) and experimentally feasible. 

6 Comments and outlook 

In this paper we have constructed different families of free particle states 
mapping harmonic oscillators states through the QAT. It would also be pos- 
sible to generalize the construction to other important families of harmonic 
oscillator states, Gaussian states, which are mixed states with the property 
that their Wigner distribution functions are Gaussian (see, for instance, 
[IT]). They include coherent states, squeezed states and thermal equilib- 
rium states. Their main property is that Gaussian states keep this property 
under evolution even in the presence of dissipation and decoherence. That 
is the reason why they have been profusely used in Quantum Optics and 
Quantum Information in the Continuous Variables formalism |18|. [19] . 

We should stress that, by using the QAT, the density matrix p' of a 
mixed harmonic oscillator state can be mapped into the density matrix p of 
a mixed free particle state: 

p = Ap'tf (41) 

The unitarity of A guaranties that p is a proper density matrix, provided 
that /jho is. Even more, it is easily checked using the evolution operator 
constructed in [TJ that: 




i.e., if pi satisfies the quantum Liouville equation for the harmonic oscillator 
Hamiltonian, then p satisfies the free particle counterpart. 
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All the properties of the density matrix p' are transferred to p, such as 
characteristic functions, quasi probability distributions, etc. In particular, 
if p' describes a Gaussian state, also p corresponds to a Gaussian state. 

A deeper study of the QAT applied to mixed states and how it can be 
used to describe dissipation and decoherence analyzing the transformation 
properties of master equations like the Lindblad one under the QAT is the 
object of work under progress and will be presented elsewhere. 

Another interesting point to study is whether the free particle states with 
many humps are physically observable and measurable, since their harmonic 
oscillator counterpart, the Hermite-Gauss and Laguerre-Gauss states are 
"visible" in Quantum Optics using lasers [5]. Let us consider, for instance, 
a two-hump wave packet <frQ' ''"' ' (x,t) in two or three dimensions with the 
humps in the transversal direction to that of the mean velocity vq. The 
separation of the two maxima of \<f>^'°' '"' | 2 (see Fig. [1]) is greater, in a 
factor 1.6, than the uncertainty in position Ax±. Therefore the two humps 
should be measurable, and in fact, if this wave packet propagates in a bubble 
or wire chamber, two parallel, divergent tracks would be observed (if times 
t « t are considered). For a three-hump wave packet, the separation 
among consecutive maxima (see Fig. [1]) is smaller than the uncertainty 
in position, although the distance between the more separated maxima is 
greater than the uncertainty in position. This, together with the fact that 
the central maximum is smaller than the external ones, suggests that only 
two, overlapping thick tracks would be observed in a bubble or wire chamber. 
A similar behavior for a larger number of humps is expected. We think that 
these points deserve further study and clarification. 

The ideas developed here could also be applied to relativistic systems, 
particularly to the free particle in de-Sitter space-time, where the ordinary 
formulation of quantum theory does not find a natural physical vacuum 
|22j . In this sense, the generalization of our approach to the relativistic 
case would provide a hierarchy of states where the first state, the relativistic 
counterpart of the Gaussian wave packet [231 Elj , plays the role of a vacuum 

m- 
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Appendix: Derivation of coherent, squeezed, num- 
ber states for the free particle 

In this Appendix we write down the configuration space wave functions for 
coherent, squeezed, number states for the harmonic oscillator, and derive 
the corresponding wave functions for the free particle using the QAT. Let 
us first introduce some notation in order to simplify the expressions. Let 
8' = 1 + ztan(wt')) and K = 1 + ie 2r tan(W)- These expressions, under 
the classical Arnold transform @, are mapped to their corresponding un- 
primmed versions, 5 = 1 + iut and S r = 1 + ie 2r 'ujt. 
Let us also define the dimensionless quantity 

(x'-XocosM') - -22- sin(wf)) 



H (e 2r sin 2 (u^) + e" 2r cos 2 (wt')) 1/2 

frnw\$'\ r( , . Po . / j» ,„„x 

= \ —r- T77T e \ x — xocos(wr) smlwi ) ) . (43) 

V ft |o,'.| raw 

The configuration space wave functions are then written as: 



V&r^.O (jfj)' (ij)^" ^^^nfoO, (44) 



where the phase 6'(x',t') is given by: 

1 



V,0 



/|2 



I* 



- W tan(wt ) g 

Vmlu V ft / 



+ sinh(2r) tan(^')<? /2 



(45) 



To compute the QAT of these states, ifi? r Jx,t), it is useful to write 
down a "dictionary", completing 

ut' ->■ tan _1 (a;t) x' — >■ t|t 

cos(wt') -> ^ sinMO -> g 

5' -> 5 5; -> 5 r (46) 

• * ^ 2 



(47) 



Using this dictionary, the adimensional quantity q' changes to: 



x ~ x o + m* 



m 



^/2Le- r |5 r | ' 
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and the phase 9(x',t') changes to: 



(./•'./') — (H.r.l) - — 12 



1/1 n mUJ '» 



+e \d r \ [—===- J— x 0Jt]q 
,\mfujj V " 



+ smh(2r)uitq 



(48) 



Therefore, the corresponding coherent, squeezed, number state for the 
free particle reads: 

^a,r)M) =^/ 2 -^ Q^J ^^)^E B ( g ). (49) 
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